Mon. Not. R. Astron. Soc. 000,[T]{To|(2008) Printed 19 March 2009 (MN KTeX style file v2.2) 



An HST/ACS View of the Inhomogeneous Outer Halo of M31* 

q< ; J. C. Richardson 1 !, A. M. N. Ferguson 1 , A. D. Mackey 1 , M. J. Irwin 2 , S. C. Chapman 2 , 
S : A. Huxor 3 , R. A. Ibata 4 , G. F. Lewis 5 , N. R. Tanvir 6 

1 Institute for Astronomy, University of Edinburgh, Royal Observatory, Blackford Hill, Edinburgh EH9 3HJ, UK 
2 Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge, CB3 OHA, UK 
3 Department of Physics, University of Bristol, Tyndall Avenue, Bristol, BS8 ITL, UK 
4 Observatoire de Strasbourg, 11, rue de V University, F-67000 Strasbourg, France 
5 Sydney Institute for Astronomy, School of Physics, A29, University of Sydney, NSW 2006, Australia 
6 Department of Physics and Astronomy; University of Leicester, LEI 7RH, UK 



Accepted ?. Received ?; in original form ? 



ABSTRACT 

We present a high precision photometric view of the stellar populations in the outer halo 
of M31, using data taken with the Hubble Space Telescope Advanced Camera for Surveys 
(HST/ACS). We analyse the field populations adjacent to 1 1 luminous globular clusters which 
sample the galactocentric radial range 18 < R < 100 kpc and reach a photometric depth of 
~ 2.5 magnitudes below the horizontal branch (nipsuw ~ 27 mag). The colour-magnitude 
diagrams (CMDs) are well populated out to ~ 60 kpc and exhibit relatively metal-rich red gi- 
ant branches, with the densest fields also showing evidence for prominent red clumps. We use 
the Dartmouth isochrones to construct metallicity distribution functions (MDFs) which con- 
firm the presence of dominant populations with < [Fe/H] >« —0.6 to —1.0 dex and consider- 
able metallicity dispersions of 0.2 to 0.3 dex (assuming a 10 Gyr population and scaled-Solar 
abundances). The average metallicity over the range 30-60 kpc is [Fe/H]= —0.8 ± 0.14 dex, 
with no evidence for a significant radial gradient. Metal-poor stars ([Fe/H]^ —1.3) typically 
account for < 10 — 20% of the population in each field, irrespective of radius. Assuming our 
fields are unbiased probes of the dominant stellar populations in these parts, we find that the 
M3 1 outer halo remains considerably more metal-rich than that of the Milky Way out to at 
least 60 kpc. 

Key words: galaxies: evolution — galaxies: formation — galaxies: halo — galaxies: individual 
(M31) — galaxies: stellar content — galaxies: structure 



1 INTRODUCTION 

Understanding the nature and origin of stars at extreme distances 
from galaxy centres is of paramount importance. Modem theo- 
ries predict that the stellar halos of massive galaxies result from 
the continu ed accretion and disru ptio n of smaller satellite sys- 
tems [from ISearle & ZinrJ Jl978h to iBullock & JohnstorJ 120051 ; 
|Pe Lucia~& Helmi 2008y. The stellar populations these systems do- 
nate to their host halos are unique in terms of their mix of ages, 
chemical compositions and kinematics. In the outer halos of galax- 
ies, whe re the mixing timesc ales are of the order of several gi- 
gayears jjohnstonet afl ll996). fossil signatures of these accretion 

* Based on observations made with the NASA/ESA Hubble Space Tele- 
scope, obtained at the Space Telescope Science Institute, which is operated 
by the Association of Universities for Research in Astronomy, Inc., under 
NASA contract NAS 5-26555. 
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events are expected to be well-preserved and detailed studies of the 
stellar populations in these parts can constrain the number, nature 
and timescales of b oth recent and ancient accretion events (e.g., 
I Johnston et afll2 008). Such information is required in order to rig- 
orously test hierarchical growth models against the more traditional 
mon olithic collapse mod els, in which halo stars form entirely in situ 
(e.g jEggen et"al]|l962h . 

The last decade has seen tremendous advances in our under- 
standing of our nearest large neighbour, M31. A particular high- 
light has been the first detailed exploration of the large-scale stellar 
distribution and content outside the main body of the galaxy. The 
Isaac Newton Telescope wide-field imaging survey of M3 1 mapped 
a ~ 100 x 100 kpc 2 region down to » 3 magnitudes below the 
red giant branch (RGB) tip and revealed copious substructure in 
the spatial distribution of stars dlbata et al] l200ll : iFerguson et al.l 
I2002L llrwin et alj|2005h . Follow-up work with the Hubble Space 
Telescope Advanced Camera for Surveys (HST/ACS) and the 
DEIMOS spectrograph on Keck has shown that almost all of this 
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inner halo substructure can be explained by a c ombination of gi 



200 



ant stream debris an d pert urbed disk material dlbata et al 
20051: iFersuson et al.ll2005l: iBrown et alJl2006t iFardal et alJl200' 
Gilbert et alJl2007l : lFaria et al.ll2007l : lRichardson et alj|2008l) . sug- 
gesting that there has only been one significant accretion event onto 
the inner halo of M3 1 over the last several gigayears. 

Attention has recently been turned to the more remote regions 
of M31's halo. Using MegaCam on the Canada-France-Hawaii 
Telescope (CFHT), efforts are underway to map the stellar dis- 
tribution out to a radius of 150 kpc around M31, with an ex ten- 
sion towards M33 dlbata et alfeoollMcConnachie et af]|2008h . In 
addition to the discovery of several new dwarf galaxies and halo 
globular clusters (GCs) dMartin et al.l 12006]: iMackev et al.l l200d 
l2007l ; lHuxor et alj2008l ; lMcConnachie et al.l2008l) . this survey has 
also uncovered new faint debris streams as well as an extremely 
extended stellar halo with a shallow radial fall-off of E^(_R) oc 
R~ 19 . Spectroscopy of individual RGB stars in these outlying re- 
gions has yielded the first evidence for a pre ssure-supported metal 
poor ([Fe/H]~ —1 . 4) population in t he halo dChapman et al .120061 : 
iKalirai et af]|2006l : iKoch et al.ll2008l) . although it remains unclear 
whether this component - as opposed to a more metal-rich popu- 
lation associated with discrete substructures - dominates the stellar 
halo mass in these parts. 

In this paper, we present a high precision view of the stel- 
lar content in the outer halo of M31 using data taken with the 
HST/ACS. We analyse the field populations adjacent to 11 lumi- 
nous GCs which lie at projected radii of 18 kpc ^ R ^ 100 kpc. 
The HST/ACS data yield accurate photometry with high complete- 
ness to ~2.5 magnitudes below the horizontal branch and thus en- 
able the construction of colour-magnitude diagrams (CMDs) and 
metallicity distribution functions (MDFs) that are a vast improve- 
ment over those obtained from ground, where Galactic foreground 
and background galaxy contamination pose significant challenges 
at even moderately bright magnitudes. 

The paper is outlined as follows. In § 2 we describe the ob- 
servations and subsequent treatment of the data while in § 3 we 
compare and contrast the detailed CMDs of the fields and derive 
their metallicity distribution functions. In § 4 we discuss our results 
in the context of theoretical predictions and current literature and 
in § 5 we summarise our findings. 



2 OBSERVATIONS AND DATA REDUCTION 

The observations were obtained with the Wide Field Channel of the 
HST/ACS as part of the Cycle 13 program GO 10394 (PI Tanvir). 
The original aim of this program was to obtain deep high resolution 
photometry of outer halo GCs in M31 that were discovered in the 
cours e of the the Isaac Newton Telescope and CFHT MegaCa m sur- 
veys dHuxor et al.ll2005l : lMartin et al.ll2006| ; lHuxor et all2008t see) 
and results f rom this analysis are reported elsewhere dMackev et al.l 
2006. 120071 : Tanvir et al. in prep.). In the present paper, we focus 
instead on the field populations imaged alongside the GCs. 

Fields were observed in the HST/ACS F606W and F814W fil- 
ters for ~ 1800s and ~ 3000s respectively, with small dithers be- 
tween various sub- exposures. Ful l details of the photometric reduc- 
tion are provided in lMackev et alj d2006t). In brief, w e used the ACS 
module of the DOLPHOT program ( Dolphin 2000) to obtain point- 
spread function (PSF) fitting photometry for each image. Next, 
quality information provided by DOLPHOT was used to clean the 
source catalogue of any non stellar sources such as background 
galaxies and blended stars. Specifically, tests showed that strong 
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Figure 1. A schematic map of M31 with the globular cluster positions 
marked as black stars. The inner ellipse represents M31's main stellar disk 
(i = 77° and R = 2° or 27 kpc) while the outer ellipse has a radius of 55 kpc 
and is flattened to b/a = 0.6. The outer circle marks a radius of 100 kpc. 
Cartoons of the most significant outer halo substructures discovered to date 
are lab elled. The inner halo substructures examined in Richardson et al. 
12008) are highlighted with unlabelled dashed ellipses. Moving from the 
northeast to the southwest, they are the NE Clump, N Spur, NGC205 Loop, 
Claw and Gl Clump. 



stellar detections satisfied the conditions —0.5 < sharpness < 0.9, 
— 1.7 < roundness < 1.7 and signal-to-noise ratio > 4. The data 
were further pruned according to the distribution of the magni- 
tude error, x 2 °f me PSF, and sharpness as a function of mag- 
nitude. Finally, magnitudes were transformed to the VEGAmag 
scale using the zero-points of ISirianni et alj d2005l) and were cor- 
re cted for foreground reddening by interpolating within the maps 
of lSchleeel et al.ldl998i) . 

Thirteen GC fields were imaged in total. Two of these sam- 
ple rich substructure i n the inner halo of M31 and have previ- 
ously been presented in iRichardson et al.1 d2008l) : we will not dis- 
cuss these fields further here. Extensive artificial star tests on those 
fields demonstrated they were > 80% complete at rriFmewfi = 
27.3 mag and m_F8iw,o = 27.0 mag. Given that they are a fac- 
tor of ~ 100 times more crowded than fields analysed in this pa- 
per, these values serve as a lower limit to our completeness. Stars 
belonging to the GCs have been excluded from our analysis by 
masking a circular region of radius 40 — 50" around each GC 
(see TableQ}. This corresponds to one tidal radius (Rt) for GC10 
and at least 2Rt for all other clusters (Tanvir et al. in prep). Ta- 
bleQ]lists the location and projected radius from the centre of M31 
for each pointing, as well as the reddening value derived from the 
ISchlegel et alj J 19981) maps. Field names follow the nomenclature 
of lMackev et al.1 d2006t [20071) . For example, GC1F is the field pop- 
ulation around the classical globular cluster GC1, and EC2F is the 
field population associated with extended globular cluster EC2. 

Figure 1 shows the locations of the GC fields around M31. 
All but one of the fields probe what we will refer to as the outer 
halo (R > 30 kpc) of M31. This boundary marks a key transition 
point in the radial surface brightness profile of the galaxy where 
the behaviour switches from a steep inner fall-off (often described 
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Table 1. Observational Information 



Field 


R.A. 


Dec 


R" 


h(B — V ) 


Rcut C 


F f9 


< [Fe/H\ > e 


17 f 




(J2000.0) 


(J2000.0) 


(kpc) 




(") 


(%) 


(dex) 


(%) 


GC7F 


00:38:49.4 


42:22:48.0 


18.2 


0.062 


40 


0.3 


-0.78 


12.0 


GC3F 


00:30:27.3 


41:36:20.4 


31.8 


0.075 


45 


2.7 


-0.82 


17.0 


GC2F 


00:29:44.9 


41:13:09.8 


33.4 


0.073 


45 


3.3 


-0.87 


21.4 


EC2F 


00:42:55.1 


43:57:28.5 


36.8 


0.087 


50 


1.0 


-0.81 


15.5 


GC8F 


00:54:25.0 


39:42:55.5 


37.1 


0.051 


45 


3.8 


-0.60 


12.6 


GC9F 


00:55:44.0 


42:46:16.1 


38.9 


0.099 


45 


1.6 


-0.67 


10.0 


GC1F 


00:26:47.8 


39:44:45.5 


46.4 


0.070 


40 


10.3 


-0.67 


10.3 


GC4F 


00:31:09.9 


37:53:59.7 


55.2 


0.058 


45 


4.7 


-1.03 


34.6 


EC4F 


00:58:15.5 


38:03:01.1 


60.0 


0.049 


40 


2.5 


-0.78 


15.5 


GC5F 


00:35:59.7 


35:41:03.6 


78.5 


0.065 


45 


44.0 






GC10F 


01:07:26.4 


35:46:49.7 


99.9 


0.055 


45 


16.8 







"Projected radial distance assumes £*M3i = 785 kpc iMcConna chie et alj2005t). 
E(B — V) values are interpolated from the reddening map of lSchlegeletalJh998h . 
c Rcut is the radial cut imposed to exclude GC stars. 

d Fj g is the percentage of foreground stars predicted in the region of the CMD used to compute the MDF (through large-area 
realisations of the Besancon model). 

e < [Fe/H] > is the average metallicity calculated from the MDF (Section 3.1). 

f F mp is the metal-poor fraction, given by relative number of stars in the MDF that have ([Fe/H]^-1.3 dex). 



with an R 1 / 4 law) to a shallower power-law form ilrwin et alj2005i 
Ibata et al. 2007; hereafter 107). Using the extended disk parame- 
ters of PA= 38.1° and i = 64.7° from Ibata et al. (2005), the 
field locations analysed here correspond to de-projected radii of 
40 — 230 kpc which makes it highly unlikely that they contain any 
substantial disk component. 



3 ANALYSIS 

3.1 Colour-Magnitude Diagrams (CMDs) 

Figure [2] shows the CMDs plotted in order of increasing projected 
galactocentric radius. The average photometric errors returned by 
DOLPHOT are indicated towards the left hand side of the pan- 
els. Theoretical isochrones of [Fe/H]=— 1.5, —0.7, —0.2 dex for 
[tv/Fe]= and age 10 Gyr have been o btained from the D art- 
mouth Stellar Evolution Database (DSED; IBotter et al J|2007h and 
are superimposed. They ha ve been shifted to the dist ance modulus 
of M31 (m - M = 24.47: lMcConnachie et alj2005l) . 

The CMDs range from being well-populated (~ 7000 stars at 
R = 18 kpc) to being very sparse 600 stars at R ^ 80 kpc). 
As a result, it is prudent to first consider the level of contamina- 
tion that might be present from the foreground Milky Way (MW) 
population. Since our fields span almost 10 degrees in Galactic 
latitude, the foreground population towards each field is expected 
to vary substantially depending on the proximity of the line-of- 
sight to t he Galactic Plane. We have used the Besancon Galac- 
tic model dRobin et alj|2003l) to obtain large area realisations (100 
x larger than the ACS field of view) of the foreground popula- 
tion towards each field which we subsequently scale down to the 
size of an ACS pointing in order to reduce the shot noise in a 
given simulation. The error function measured in each field was 
fed into the model and the results were de-reddened in the same 
way as the data. We find the number of foreground stars predicted 
to lie in the region of the CMD used to derive the MDFs ("^ 13.2b is 
less than or equal to 10% in all but the two outermost fields (see 
Ff g in Table 1). Furthermore, we find that the colour and magni- 
tude distribution of foreground stars is unlike that of our observed 



CMDs. The bottom right-hand panel of Figure 2 shows the CMD 
of a large area realisation of the Besancon model towards a rep- 
resentative field (GC2F) for increased resolution. This CMD con- 
tains two orders of magnitude more stars than are predicted to fall 
within a single ACS field and is intended as a comparison only. The 
CMD of the foreground clearly does not reproduce any of the dis- 
tinct features observed in the GC fields, such as the red clump at 
m F(i o6w,o - m F8 i4w,o = 0.8 mag, m F 8i4w,o = 24.0 mag, nor 
does it follow the curvature of the RGBs. We can therefore be con- 
fident that foreground contamination is not a major issue for most 
of our fields. 

Figure [2] shows prominent RGB sequences in all fields out to 
R = 60 kpc with the bulk of the stars well-bounded by the [Fe/H] 
= -1.5 dex and -0.2 dex isochrones. The mean RGB loci are consis- 
tent with a moderately high metallicity of [Fe/H]~ —0.7 although 
one field, GC4F, appears more metal poor with most stars lying 
between the [Fe/H]= —1.5 and —0.7 dex isochrones. The widths 
of the RGBs are considerably larger than the photometric error bars 
for magnitudes brighter than m,F8i4W,o ~ 26.0 mag suggesting an 
intrinsic spread in metallicity. It is worth considering whether line- 
of-sight effects could artificially create this width. In overlaying the 
isochrones, we have assumed that the stars in each field all lie at the 
systemic distance of M31, whereas in reality they will span a range 
of line-of-sight distances. We checked the effect a distance spread 
would have on the CMDs by plotting isochrones with varying dis- 
tance moduli and measuring the resulting colour width. In the ex- 
treme case of shifting the —0.7 dex isochrone to lie ±60 kpc in 
front (behind) M31, this induced a colour broadening of 0.02 mag 
at m_F8i4W,o = 25.0 mag and 0.05 mag at rriFsi4w,o ~ 22.0 mag. 
This accounts for less than 10% and 15% of the width of the RGB 
at those magnitudes respectively and implies the RGB spreads are 
for the most part intrinsic. 

A red clump (RC) due to core helium burning stars is detected 
in most fields with a well-populated RGB (GC7F, GC3F, GC2F, 
EC2F, GC9F, EC4F) and indicates the presence of an inter mediate 
age and/or metal-rich population dGirardi & Salarisll200ll) . There 
is a tentative detection of an old, metal-poor extended horizon- 
tal branch (HB) in a few fields (e.g. GC7F, GC3F, EC2F, GC9F) 
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Figure 2. The extinction-corrected CMDs of the outer halo fields. Listed below the fi eld name are the projected radial distance and the total number of sources 
on the CMD. Isochrones from the Dartmouth Stellar Evolution Database (DSED; lDotter et alj|2007l) with age 10 Gyr and [Fe/H ]=-1.5 (green), [Fe/H]=-0. 7 
(blue) and [Fe/H]=-0.2 (red) have been plotted for comparison. They have been shifted to the distance of M31 (m-M)=24.47 ( McConnachie et al. 2005). 
Typical photometric errors (grey crosses) are plotted. GC5F and GC10F do not show convincing RGBs but show a clustering of sources at faint magnitudes. 
Other fields have CMDs whose density at the faint end scales with the density of the red clump suggesting that some fractio n of these sources are genuine stars. 
The last panel shows the MW foreground (FG) contamination predicted by the Besancon Galactic model jRobin et al.l2003l) over an area 100 times larger than 
our ACS fields (centred on GC2F, a representative field). This large area realisation is included to produce a smooth representation of the foreground CMD 
morphology to compare with the GC fields. 
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which can be traced from mFeoewfi — wiF8i4w,o = 0.50 mag, 
mF8i4w,o = 24.24 mag to rn F 606w,o - m F si4.w,o = 0.00 mag, 
tof814w,o = 25.75 mag. This feature may also be present in GC8F 
and EC4F but there are too few stars detected to be certain. Indeed, 
the sparseness of most of the CMDs makes it impossible to do com- 
parative ana lyses of the RC/HB mor phologies in different fields, as 
was done in Richardso netai] J2008I) . 

Figuref2]also highlights the presence of star count (i.e. surface 
brightness) and stellar population inhomogeneities throughout the 
halo. For example, EC2F, GC8F and GC9F lie in the projected ra- 
dial range of ~ 37 — 39 kpc yet their total star counts fluctuate by 
a factor of three (2437 ± 49, 862 ± 29 and 1445 ± 38 respectively) 
indicating genuine variations larger than the Poisson uncertainties. 
Furthermore, EC4F and GC4F lie in the projected radial interval 
55 — 60 kpc, however EC4F exhibits a wider RGB than GC4F. 
The former also contains twice as many stars (1138 ± 34) as GC1F 
(666±26), a field which lies 15 kpc further in. Since the fields have 
a uniform depth and completeness these patterns indicate the that 
they do not probe a smooth well-mixed stellar component of M3 1 . 

Although we detect stars in fields GC5F (R = 79 kpc) 
and GC10F (R = 100 kpc), the CMDs are too sparsely pop- 
ulated to reliably constrain their stellar populations. Most of the 
sources in these CMDs cluster around the base of the RGB se- 
quences and while at least some of them will be genuine RGB stars 
(note that in the more populated fields, the density of stars with 
TiF8i4W,o < 26.0 scales with the RGB and RC), there will also 
be a high contribution from unresolved background galaxies and 
foreground stars. Without knowledge of the true background popu- 
lations along these sight-lines, it is impossible to disentangle these 
different components at faint magnitudes (mF8i4w,o > 25.5 mag). 

3.2 Metallicity Distribution Functions (MDFs) 

In order to further characterise the stellar populations in each field, 
we determi ne the MDFs usin g a set of theoretical isochrones from 
the DSED l lDotter et al1l2007l) with age 10 Gyr and a Solar-scaled 
abundance mixture. Figure [3] shows the RGB sequences of the se- 
lected isochrones shifted to the distance of M3 1 and overlaid onto 
the photometry of EC2F. We have selected isochrones which form 
a finely-spaced grid ranging from [Fe/H]= —2.5 dex to +0.3 
dex in [Fe/H]= 0.1 dex steps. With such sampling, the metal- 
licities of RGB stars can simply be derived by comparing their 
TO-F606W.0 — rriF8i4w,o, TiF8i4w,o values to the isochrone grid 
and assigning the metallicity of the closest isochrone. 

In order to minimise contamination, the MDFs are computed 
from RGB stars selected within the range 21.0 < raF8i4W,o < 
25.25, mF606w,o — rriFsawfi < 1.7 and sandwiched between 
the -2.5 and +0.3 dex isochrones. By selecting stars brighter than 
TiF8i4w,o = 25.25 mag, we limit the influence of unresolved 
background sources and ensure a minimum completeness rate of « 
95%. This lower magnitude limit marks the 'knee' in the m F8i4W,o 
luminosity function of the co-added fields beyond which the source 
count rises very rapidly. The upper cut-off of mF8i4iv,o = 21.0 
mag minimises the inclusion of bright asymptotic giant branch 
(AGB) stars and the colour cut (mF606W,o — rriFsuw.o < 1-7) 
avoids areas where the density of stars on the CMD is very small 
and therefore liable to increased influence from the MW fore- 
ground. An additional area (23.9 < mF&uwfi < 24.6) was cut 
to exclude the RC region. Although RGB stars fainter than the RC 
have less metallicity sensitivity than their brighter counterparts, 
they represent the bulk of the RGB populations in these sparse 
fields and are essential for providing a sufficient sample of stars. 
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Figure 3. The CMD of EC2F (grey dots) on top of w hich is shown the 
theoretical RGB sequences of DSED iDotteret all 2007) used to compute 
the metallicity distributions. The isochrones have a Solar-scaled abundance 
mixture, age 10 Gyr and span the range [Fe/H] = —2.5 to +0.3 in 0.1 dex 
steps (left to right). They have been shifted to (m-M) = 24.47. 

We have verified that our results do not change significantly if we 
only use stars exclusively above or below the RC. 

We present the resultant MDFs of the fields in Figure|4](solid 
black lines). The average metallicity and standard deviation of each 
field are quoted in the panels and marked on the MDF by a black 
dashed line (see also Table 1). Out to R = 60 kpc, the MDFs 
are generally characterised by broad distributions which peak at 
intermediate metallicities < [Fe/H] >~ —0.6 to — 1.0 dex, con- 
firming inferences from visual inspection of the CMDs alone. We 
reiterate that some of the metallicity spread might be due to pho- 
tometric errors and line-of-sight depth but these will not domi- 
nate (see the top panel of Figure 4). Averaging over GC fields in 
the 30-60 kpc range, the mean stellar metallicity is found to be 
< [Fe/H] > = —0.78 ± 0.14. As previously mentioned, there is no 
significant RGB detected in GC5F and GC10F and hence we are 
unable to compute individual MDFs for these fields. However, they 
provide a useful probe of the 'contamination' in our MDF selection 
box from non-M31 sources, namely unresolved background galax- 
ies and foreground dwarf stars. The bottom left-hand plot shows 
the MDF derived for the average of the GC5F and GC10F fields. 

Also shown in Figure|4]are solid grey lines representing MDFs 
for the MW foreground t owards each field as predicted by the 
Besancon Galactic model l lRobin et alj|2003h . These are computed 
by taking the simulations described previously and computing an 
MDF in an identical manner to our GC fields. These are obviously 
not true representations of the metallicities of MW foreground stars 
since this population is dominated by dwarfs and not RGB stars on 
which our analysis is based. Nevertheless, this procedure provides 
additional insight into the degree to which the foreground might 
be biasing our results. The average metallicity of the foreground is 
listed in grey and highlighted with a grey dashed line. The fore- 
ground MDF peaks at high metallicities ([Fe/H]~ —0.4 dex) and 
has a long metal poor tail, a shape which is rather different from 
that observed in any of the fields out to 60 kpc. Thus, foreground 
stars alone cannot explain the observed form of the MDFs with 
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Figure 4. The MDFs for our outer halo fields assuming an age of 10 Gyr and a Solar-scaled abundance distribution. Reported in each panel are the name 
and projected radial distance of the field plus the number of stars used to compute the MDF. The average metallicity < [Fe/H] > and its associated standard 
deviation, are listed below the field name and marked on the plot with a black dashed line. Gray lines and text indicate results for the foreground MW 
population. The horizontal bars in the top left plot represent the typical error in the metallicity derived from the colour uncertainties at rrii?8l4VK0 = 25.0 
mag. 
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Table 2. The effect of varying age and [ct/Fe] on the derived mean metal- 
licity. 



Age 


[a/Fe] 


< [Fe/H] > shift 


(Gyr) 


(dex) 


(dex) 


12 


0.0 


-0.01 


10 


0.0 


0.00 


10 


0.2 


-0.04 


8 


0.0 


+0.02 


6 


0.0 


+0.05 



their strong peaks at intermediate metallicities (although the fore- 
ground could explain a large fraction of the Solar metallicity stars 
observed). We examined the impact of subtracting the MW fore- 
ground MDFs from each field and found it to be minimal. The 
average metallicities decreased by only ~ 0.02 dex in all fields 
except for GC1F, which became 0.08 dex more metal-poor, while 
the metallicity dispersions remained essentially unchanged. 

The shape of the average MDF of the GC5F+GC10F fields 
and the Bescancon model prediction are generally in good agree- 
ment which suggests a large fraction of the contamination in the 
MDFs is due to foreground stars. To test this further we performed a 
two-sided Kolmogorov-Smirnov (KS) test and found a high proba- 
bility (P = 0.95) that the MDFs of GC5F+GC 10F and the Besancon 
foreground are drawn from the same parent distribution. Note that 
this is true only for bright stars which were incorporated into 
the MDFs. At magnitudes fainter than tufsiaw = 25.25 GC5F 
and GC10F suffer additional contamination from unresolved back- 
ground galaxies. The surplus of stars at [Fe/H]< —0.8 dex here 
may be a mixture of bright unresolved background sources and true 
metal-poor M31 halo stars, or else it may highlight a slight normal- 
isation problem with the Galactic model. 

Many of the MDFs contain stars with metallicities reach- 
ing as low as —2.5 dex (this being the lower limit considered 
in our MDF derivation). Interestingly, the fraction of metal-poor 
([Fe/H] ^ —1.3 dex) stars (F mp ) is relatively invariant from field 
to field at ~ 10 - 20% although one field, GC4F at 55 kpc, has a 
slightly larger metal-poor fraction (see Table 1). Since it is possi- 
ble that some of this population could be due to AGB stars on the 
blue edge of the RGB, this number provides an upper limit to the 
fraction of metal-poor halo stars sampled by our fields. 

In deriving the MDFs, we have made the simplifying assump- 
tion that the stellar populations are uniformly 10 Gyr old. The ben- 
efit of choosing an old age to compute an MDF is that the posi- 
tion of an old star within the RGB is mostly determined by its 
metallicity and largely insensitive to its age. In reality, our fields 
will contain com posite stellar po pulations with varying age and 
metallicity (e.g., Bro wn et al. 2008). We have examined the ef- 
fect of a straight scaling of the age and alpha-enhancement of the 
isochrones on the derived MDF (Table 2). Specifically, we recon- 
structed the MDF of each field exactly as before but with a new 
set of isochrones as described in Table 2. The mean shift in aver- 
age metallicity (< [Fe/H] > shift) was calculated from the resulting 
MDFs of all fields. When isochrones of 6 Gyr are used, there is only 
a modest +0.05 dex increase in the average metallicity. On the other 
hand, increasing [a/Fe] to +0.2 has the effect of shifting the mean 
metallicity by -0.04 dex. While the assumption of a fixed age and 
alpha-enhancement will undoubtedly introduce uncertainties into 
our results, we do not believe that these will be significant. 



4 DISCUSSION 

Our analysis of the stellar populations in the vicinity of M3 1 outer 
halo GCs indicates a field-to-field variation in the mean metallicity 
of < [Fe/H] >~ —0.6 to —1.0 dex with a typical dispersion of 0.2- 
0.3 dex. These relatively high metallicities are supported by both 
CMD morphologies (e.g. mean RGB colours and the presence of 
prominent RCs) as well as the MDFs. The average metallicity over 
the range 30-60 kpc is [Fe/H]= — 0.8±0.14 dex, with no evidence 
for a significant radial gradient. 

Our findings for the M31 outer halo are in good agree- 
ment with expectations from models in which stellar halos are 
built predominantly from satellite accretion jFont et alj2003.l2008l : 
I5e Lucia & Helmj|2008h . These models predict stellar halos to be 
highly structured in their outskirts with discrete substructures be- 
coming increasingly dominant at large radius, an effect mostly 
driven by the long mixing times in these parts (e.g.. |johnston et all 
2008). While the pencil beam observations discussed here cannot 
test how the fraction of substructure increases with radius, the field- 
to-field variations we observe in RGB star counts provide good sup- 
port for the notion that the outer regions are poorly mixed in terms 
of stellar density. 



4.1 Connections to Known Substructures 

Even though our fields were selected purely on the basis of their 
encompassing a luminous outer halo GC, Figure 1 reveals that a 
good fraction of them lie on, or near, discrete substructures iden- 
tified in the CFHT MegaCam survey. This is intriguing but per- 
haps not completely surprising given that a significant fraction of 
the halo luminosity out to 60 kpc resides in such structures, which 
have a large covering area on the sky (107). GC7F, GC3F and 
GC2F are situated in reasonable proximity to the diffuse West- 
ern shelf and GC9F is near to the North-Eastern Shelf. Both of 
these substructures are most likely comprised of debris stripped 
from the metal-rich g iant stream progenitor {Fardal et alj 1 20071 : 
iRichardson et alj 2008). GC8F is located in close proximity to the 
minor axis transverse Stream D, while EC4F probes Stream C. 
IChapman et alj {2008) have recently found a bimodal spectroscopic 
metallicity distribution for Stream C with peaks at [Fe/H]=— 0.7 
and [Fe/H]=— 1.3 dex. Our MDF for EC4F also supports a bimodal 
distribution with peaks at [Fe/H]=— 0.85 and [Fe/H]=— 1.45 dex, in 
good agreement. iFardal et alj 120081) have shown it is plausible that 
many of these outer halo features originate from the same event, 
namely the accretion of the (metal-rich) giant stream progenitor. 
Alternatively, even if the outer halo debris originates from several 
different accretion events, the spatial coherence and moderately 
high surface brightness of the streams imply th at the progenitor s 
were of intermediate mass and hence metallicity {Font et al. 2008). 

The most metal-poor population identified in our study is 
GC4F. This field probes the 'major axis diffuse structure' dis- 
covered by 107, who find a metallicity of [Fe/H]=— 1.4 dex from 
the colour of the RGB. Our MDF shows a dominant peak at 
[Fe/H]~ —1.2 dex in this structure, which is in reasonable agree- 
ment. GC4 has been recently sh own to have a population of extra- 
tidal stars dFederici et al J [20071 : Tanvir et al. in prep.). In deriv- 
ing the GC4F MDF we have taken measures to avoid contami- 
nation by this population by masking out the cluster to twice its 
tidal radius (Rt = 22.4"; Tanvir et al. in prep.) but it is worth- 
while considering whether this has been sufficient. In Figure 5 we 
compare the MDFs of the field and cluster populations where the 
latter are defined as stars lying within 3 half-light radii (where 
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GC4 field -1.03 ± 0.33 55 stors 

GC4 cluster! -1 .79 ± 0.06 17 stars 




Figure 5. The relationship between the globular cluster GC4 and the as- 
sociated field population. The solid black line represents the MDF of the 
field population (with the GC masked out to two Rt ) while the grey dashed 
line is the globular cluster MDF measured up to three half-light radii (Rh)- 
The average metallicity calculated over the whole metallicity range is listed 
beside the appropriate label and highlighted on the plot with a vertical line. 



Rh = 1.3"; Tanvir et al. in prep.). It can be clearly seen that the 
cluster and field star MDFs are very different. The cluster MDF 
shows a dominant peak at [Fe/H]~ —1.85 dex and an average 
metallicity of < [Fe/H] >= —1.79 ± 0.06 dex. This is in excel- 
lent agreement with the independently derived mean metallicities 
of [Fe/H] =—1.8 to —2.1 dex rep orted in the literature dGalleti et al.l 
l2005tl2006l : lMackev et af]|2007h . The more enriched field popula- 
tion (< [Fe/H] >= —1.03 ± 0.33 dex) therefore cannot be due to 
extra-tidal cluster stars, though some of the metal-poor tail could 
be. This lends credence to the major axis structure being a gen- 
uinely more metal-poor feature that is distinct in origin from the 
giant stream debris. 



4.2 Radial Metallicity Trends 

In Figure 6 we present the radial metallicity distribution of M31's 
outer halo as derived in this study and compare it to results 
in the literatur e . We i nclude here t h e spec t roscopic studies of 
IChapman et all J2008I). IKalirai et all J200d) . iKoch et al.l J2008h 
and Chapman et al THooi) " as well as the photometric studies by 
iBrown et alj d2006ll2008h . It is important to note that all the spec- 
troscopic studies invoke some degree of kinematic selection to iso- 
late the sample of halo stars which will somewhat complicate the 
comparison of their results to those from purely photom etric work. 
Additionally, the metallicities from lBrown et al.1 l200ej 120081) rep- 
resent average values derived from star formation history fits to the 
main sequence tum-off in ultra-deep HST CMDs, a method rather 
different from that which we have employed for our photometric 
metallicity derivation. 

Figure 6 shows that the metallicities we infer in our fields 
agree well with most other studies which have probed a similar ra- 
dial range although there is considerable dispersion at all radii. This 
dispersion may be due to different methods of sample selection 
and metallicity determination or to intrinsic metallicity variations 
in stellar populations throughout the halo. In support of the last 
point we highlight the ~1 dex spread in metallicity at 7? = 60 kpc 
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Figure 6. The mean metallicity and metallicity dispersion of the outer 
halo as a function of radius. Our results (red filled circles) are compared 
to all relevant studies in the literature. Gray squ ares represent the pho - 
tometric [Fe/H] of stellar streams measured by Chapman et al. (2008). 
Blue diamonds are the mean photometric metallicities of individual fields 
from [ K alirai et aT while the spectroscopic metallicities derived by 

IKoch et a?T ^^ ^) are shown as black asterisks. The photometric metallic- 
ities of the lBrown et~aj] <200d Exjl) fields are given by open circles. The 
green dashed ( dotted) lines represent the spectroscopic metallicity (disper- 
sion) traced bv lChapman et al.l 12006) between 10 < 7?(kpc) < 70. 



wh ere two separate strea ms (one metal-rich and one metal-poor - 
sec IChapman et ail 2008) are probed. 

Despite the general concordance between our results and 
those in the literature, our measurements are systematically 
more metal-rich than the sp ectroscopic metalli cities measured by 
IChapman etafl J2006h and I Koch et alj d2008l) . The discrepancy 
with Chapman et al. (2006) may be due to the different sample se- 
lection; in their analysis, they windowed out any stellar component 
that rotates within 100 km s _1 of the thin dis k over the radial ran ge 
10 < R(kpc) < 70. The discrepancy with lKoch etakl d2008l) is 
more p uzzling. These autho rs have re-analysed much of the same 
data as IKalirai et al.l J20061) using a similar halo sample selection 
method, but a different method of metallicity determination, and 
find a strong metallicity gradient which reaches [Fe/H]= —1.8 dex 
at R ~ 60 kpc. It is hard to reconcile such a metal-poor population 
with the CMD morphology and MDF form that we have observed 
in our outer halo GC fields. 

Considered together, the results plotted in Figure 6 suggest 
that much of the inner 60 kpc of M31's stellar halo is dominated 
by relatively metal-rich populations (—0.4 < [Fe/H] < —1.1). 
This situation appea rs vastly different from that in the Milky Way. 
ICarollo et alj d2007t) have recently derived metallicities for 1,235 
blue horizontal branch stars in the Galactic halo lying in the radial 
range 5-80 kpc and find the population gets progressively more 
metal poor with radius, with the dominant population at R> 30 kpc 
having < [Fe/H] >^ -2. While the halos of the Milky Way and 
M3 1 have been known for some time to differ in their metallicities 
at smaller radii, this is the first clear demonstration that these differ- 
ences remain out to significant radii. Unfortunately, few constraints 
exist on the detailed nature of the M3 1 stellar populations beyond 
80 kpc where very metal-poor populations were hinted at by the 
ground-based analysis of 107. It may be that only in these parts of 
M31 a metal poor stellar halo prevails. 

Finally, it is worth considering the possibility that by exam- 
ining the fields surrounding luminous GCs at large radii we have 
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introduced a selection bias. The outer halo GCs targeted in the cur- 
rent study comprise 1 1 ou t of the 40 (~ 25% ) we have identified 
in this particular quadrant fauxor et al.ll2008h . These objects were 
chosen for high resolution imagery with HST based based on the 
fact that they were known at the time of proposal writing and span 
a broad range of projected galactocentric radii; there are no signs 
that they are unusual compared to other M31 halo GCs. Within the 
Milky Way, at least eight GCs have been possibly associated with 
the disrupting S gr dwarf galaxy, accou nting for ~ 20% of all GCs 
at R > 10 kpc (Bell azzini et al.ll2003l) . This raises the question as 
to whether GCs at large galacto-centric radii are preferentially as- 
sociated with halo tidal stream debris, an intriguing idea which we 
will explore more in future work. 

5 SUMMARY 

We have presented high precision photometry of the stellar pop- 
ulations adjacent to 11 luminous outer halo GCs in M31 using 
HST/ACS data. Our fields sample the extended halo in the radial 
range 18 < R < 100 kpc and reach a photometric depth of ~ 2.5 
magnitudes below the horizontal branch (mFsuwfi ~ 27 mag). 
Eight pointings lie at R > 35 kpc and provide the first detailed 
view of the properties of halo stars in these remote parts. 

The colour-magnitude diagrams (CMDs) are well populated 
out to ~ 60 kpc and exhibit relatively metal-rich red giant 
branches, with the densest fields also showing evidence for promi- 
nent red clumps. Using the Dartmouth isochrones, we have con- 
structed MDFs which confirm the presence of dominant popula- 
tions with < [Fe/H] >~ —0.6 to —1.0 dex and considerable in- 
ternal metallicity dispersions of 0.2 to 0.3 dex (assuming a 10 Gyr 
population and scaled-Solar abundances). The average metallicity 
over the range 30-60 kpc is [Fe/H]= —0.8 ± 0.14 dex, with no ev- 
idence for a significant radial gradient. Although some metal-poor 
stars ([Fe/H]^ —1.3) are present in each field, they amount to typ- 
ically < 10 — 20% of the total population, irrespective of radius. 

Significant inhomogeneities are present in the fields in terms 
of both random fluctuations in star counts and field-to-field varia- 
tions in metallicity supporting the notion that the stellar populations 
are not well-mixed beyond 30 kpc, in agreement with results from 
wide-field surveys and predictions of models in which stellar halos 
assemble predominantly from satellite accretion. Our photometric 
metallicity measurements are in good agreement with most, but not 
all, spectroscopic studies of M31's outer halo. Assuming that our 
GC fields are unbiased probes of the dominant stellar populations 
in these parts, we find that the M3 1 outer halo remains considerably 
more metal-rich than that of the Milky Way out to at least 60 kpc. 
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